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Abstract 

Few-body systems with resonant short-range interactions display universal proper- 
ties that do not depend on the details of their structure or their interactions at 
short distances. In the three-body system, these properties include the existence of 
a geometric spectrum of three-body Efimov states and a discrete scaling symmetry. 
Similar universal properties appear in 4-body and possibly higher-body systems as 
well. We set up an effective theory for few-body systems in a harmonic trap and 
study the modification of universal physics for 3- and 4-particle systems in external 
confinement. In particular, we focus on systems where the Efimov effect can occur 
and investigate the dependence of the 4-body spectrum on the experimental tuning 
parameters. 

Resume 

Physique universelle de quelques corps dans un piege harmonique 

Les systemes a petit nombre de corps avec des interactions de courte portee reso- 
nantes possedent des proprietes universelles qui ne dependent pas des details de leur 
structure ou de leurs interactions aux petites distances. Dans des systemes a trois 
corps, ces proprietes incluent l'existence d'un spectre geometrique d'etats d'Efimov 
a trois corps et une symetrie d'echelle discrete. Des proprietes universelles similaires 
apparaissent dans des systemes a quatre corps et peut-etre aussi dans des systemes 
a plus grand nombre de corps. Nous construisons une theorie effective pour des 
systemes a quelques corps dans un piege harmonique et etudions la modification de 
la physique universelle des systemes a trois et quatre corps dans un confinement 
exterieur. En particulier, nous nous concentrons sur les systemes ou l'effet d'Efi- 
mov peut apparaitre et nous examinons la dependance du spectre a quatre corps en 
fonctions des parametres ajustables experiment alement. 
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1 Introduction 

Few-body systems close to the unitary limit show interesting universal prop- 
erties. In the 2-body system, the scattering amplitude saturates the unitarity 
bound from conservation of probability. Such systems are characterized by a 
scattering length a that is much larger than the typical range of the interac- 
tion ro and a is the only relevant length scale at low energies. If a is positive, 
two particles of mass m form a shallow dimer with energy E 2 ~ —h 2 /(ma 2 ) , 
independent of the mechanism responsible for the large scattering length. Ex- 
amples for such shallow dimer states are the deuteron in nuclear physics, the 
4 He dimer in atomic physics, and possibly the new charmonium state X(3872) 
in particle physics [T|2"] . 

In the 3-body system, an additional length scale 1/ac* is generated by the 
Efimov effect [3] • If at least two of the three pairs of particles have a large 
scattering length \a\ ^> ro, the Efimov effect can occur. In the limit l/a — > 0, 
there are infinitely many 3-body bound states with an accumulation point 
at the 3-body scattering threshold. These Efimov states or trimers have a 
geometric spectrum [3J: 

= -{e-^ /S0 ) n - n *h 2 Kl/m, (1) 

which is specified by the binding momentum of the Efimov trimer labeled 
by n*. This spectrum is a consequence of a discrete scaling symmetry with 
discrete scaling factor e 71 / 50 . In the case of identical bosons, Sq ~ 1.00624 and 
the discrete scaling factor is e 7 ^ 50 ~ 22.7. The scaling symmetry becomes also 
manifest in the log-periodic dependence of scattering observables on the scat- 
tering length a [I]. The consequences of discrete scale invariance and "Efimov 
physics" can be calculated in an effective field theory for short-range interac- 
tions, where the Efimov effect appears as a consequence of a renormalization 
group limit cycle [5]. 

While the Efimov effect was established theoretically already in 1970, the first 
experimental evidence for an Efimov trimer in a trapped gas of ultracold Cs 
atoms was provided only recently by its signature in the 3-body recombination 
rate [6]. It could be unravelled by varying the scattering length a over several 
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orders of magnitude using a Feshbach resonance. Since this pioneering experi- 
ment, there was significant experimental progress in observing Efimov physics 
in ultracold quantum gases. More recently, evidence for Efimov trimers in 3- 
body recombination was also obtained in a balanced mixture of atoms in three 
different hyperfine states of 6 Li [7|8] . in a mixture of Potassium and Rubidium 
atoms [9], and in an ultracold gas of 7 Li atoms jlOJ . In another experiment 
with Potassium atoms [llj, two bound trimers were observed. Efimov states 
can also be observed as resonances in atom-dimer scattering. Such resonances 
have been seen in atom-dimer mixtures of Cs atoms [T2] and of 6 Li atoms 
[T3"fT4"] . The first direct observation of Efimov trimers of 6 Li atoms created 
by radio frequency association was recently reported by the Heidelberg group 

One of the most exciting recent developments in universal few-body physics 
involves universal tetramer states. There is a pair of universal tetramer states 
associated with every Efimov trimer [16] . The resonant enhancement of 4-body 
recombination provides a signature for these tetramers similar to the trimers 
|17j . Loss resonances from both tetramers were subsequently observed in an 
ultracold gas of 133 Cs atoms [18]. In 7 Li atoms, even two sets of tetramers 
that are close to the corresponding Efimov trimers could be observed [T9"] . 
Moreover, there is some theoretical evidence of even higher-body universal 
states [20]. 

These experiments were carried out in a regime where the influence of the 
trap on the system could be neglected. However, the trap also offers new pos- 
sibilities to modify the properties of few-body systems. In particular, a tight 
harmonic confinement changes the spectrum and can lead to interesting new 
phenomena. The 2-body problem in an isotropic harmonic trap was solved 
analytically by Busch et al. [21]. In the unitary limit of infinite scattering 
length, Werner and Castin calculated the 3-body spectrum for two-component 
fermions and bosons [22] • In the case of bosons, they found that only a subset 
of states displays the discrete scaling symmetry. There are also a number of 
numerical studies of few-body systems with three and more particles in har- 
monic confinement. Most of this work, however, has focused on the problem of 
two component fermions in or near the unitary limit [52124 25 26 27 28 2SJ50] • 

In this work, we investigate universal few-body physics in harmonic confine- 
ment for three and four particles. Our strategy follows Stetcu et al. [31] , where 
an effective theory for short-range forces in the framework of the no-core shell 
model was formulated. The effective interactions were defined within a finite 
model space with a cutoff N on the basis functions. This strategy was also 
applied to atomic systems of three and four spin- 1/2 fermions in a trap [26]. In 
[28|29] , they found that using different cutoffs for systems with different num- 
ber of particles leads to improved convergence of perturbative higher-order 
corrections. This is because the spectator particles in a many-body system 
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can carry some of the excitation energy, leaving less available to an inter- 
acting pair. For sufficiently large cutoffs, however, the convergence problems 
disappear. Here, we stay at leading order and calculate the spectra of 2-, 3- 
and 4-particle systems in a harmonic trap. We focus on systems displaying the 
Effmov effect and perform a systematic study of universal bound state prop- 
erties. In particular, we investigate how the bound state spectrum depends on 
the 2- and 3-body energies required as input. Some preliminary results of our 
study were already reported in [32] . 

The paper is organized as follows. In the next section, we give a brief overview 
of our effective theory for few-body systems with large scattering lengths in 
an external confinement. In Sec. [3j we present our results for 3-body systems 
in the unitary limit and compare with the analytical solution of Werner and 
Castin [22]. Moreover, we calculate the energies of the first two L p = + states 
for a 3-body system of 6 Li in the lowest two spin states. In Sec. HI we present 
our results for 3- and 4-body systems away from the unitary limit in the form 
of an extended Effmov plot. We consider systems of two identical fermions 
with two other distinguishable particles as well as systems of identical bosons. 
The paper then ends with a summary and conclusions. 



2 Effective theories for trapped systems 



In this section, we formulate our effective theory for particles with large S- 
wave scattering length in a harmonic trap. Following Ref. |3f] , we identify the 
ultraviolet cutoff of the effective theory with the model space cutoff on the 
harmonic oscillator basis functions. Since the trap only modifies the infrared 
properties of the system, the renormalization in the ultraviolet remains the 
same and we require both a 2- and a 3-body interaction at leading order in 
the large scattering length (5]. Higher order corrections will be left for a future 
publication. For a discussion of higher order corrections for fermions with two 
spin states, see Refs. 



For A bodies with equal masses in an isotropic harmonic oscillator potential 
(HOP), the energy spectrum is determined by the Hamiltonian 



i=l 



Pi 



2m 



H — mm' 
2 



A A 
i<j i<j<k 



(2) 



where u> is the trapping frequency and V# and Wijk are 2- and 3-particle 
contact interactions between bodies i,j, and k. Since the interactions only 
depend on relative coordinates, Jacobi coordinates s*j are introduced and the 
dynamics of the centre of mass separate. As in free space, the ultraviolet 
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behavior of this Hamiltonian has to be regularised. Thus the Hilbert space is 
restricted by delimiting the basis functions with a regulator N. For A bodies 
in the HOP it is convenient to use the tensor product of harmonic oscillator 
functions (HOF) </» ni / <TOj (sA as a basis. Then the model space for a given N 
is the linear hull of HOF with the requirement J2t=i (^ n i + h) < N [31J . This 
means, the basis consists of states with an unperturbed eigenenergy less or 
equal to J2t=i(^ n i + h + 3/2)hu. In this sense, the regulator iV corresponds 
to an energy cutoff. We identify this cutoff with the ultraviolet cutoff of our 
effective theory. Since the model space is finite, the Hamilton matrix is finite, 
too. Thus the energy eigenvalues in the model space can be calculated by 
a numerical diagonalisation of the Hamilton matrix. Note that in free space 
one usually uses a momentum cutoff A to regularize the resulting integral 
equations [T|5] . For large values of N, A is proportional to y/~N. 



2.1 Two bodies 



In the 2-body sector the elements of the Hamilton matrix are 

[nlm HfJ nt n'l'm') = ftu ^2n + I + ^ 6 n , n >6 ltl >6 m>m i + —rfnfn'Vi',o (3) 



with /„ = J ( " !2 n ■ The running coupling constant v(N) is renormalised 
by the requirement that the ground state energy in the model reproduces a 
given value. The separability of the interaction allows one to find the relation 
between the ground state energy = tkj, the cutoff parameter N and 
the running coupling constant v(N): 

3 — 2 

__ 7l2hu ' \f n \ 



The energy spectrum of the 2-body system in an isotropic harmonic trap is 
known exactly [2TJ. In Fig. (TJ we show the energy spectrum as a function of 
the scattering length a in free space measured in units of the oscillator length 
(3 = Jh/ (mio). In the unitary limit, the ground state energy approaches ftw/2. 

We compare our numerical results to the exact solution in order to test the 
accuracy of our method. In Fig. |2j we show the results for the calculated 
spectra for = —1 and cutoff parameters N < 600. The exact results [21] 
are added at — > oo. Already for finite values of the cutoff A^ with Nftw large 
compared to a calculated energy, we observe good agreement with the exact 
solution. The exact values can be reproduced numerically by extrapolating 
the results in the model spaces. We note, however, that this extrapolation is 
not strictly necessary if we interpret Eq. (j2J) as an effective Hamiltonian that 
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Figure 1. Energy spectrum for two particles in a harmonic trap as a function of the 
scattering length a in units of the oscillator parameter j3 according to |21j . 
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Figure 2. Energy spectrum for two particles in a harmonic trap as a function of the 
cutoff N and exact results from [21 J at N — > oo for Eq = —Hu). The dashed line 
indicates the value of the cutoff ./V at which ./V ~ 5E/(Hoj). 



is only accurate up to errors of order yE/(Nhui). The dashed line indicates 
the value of the cutoff iV at which iV »s 5E/(hu>). If this value of N is reached, 
the effects from the finite cutoff are small and the extrapolation is usually not 
required. If more bodies are considered, however, it is computationally more 
expensive to go to large values of iV and extrapolation becomes necessary. In 
order to test the accuracy of our method, we have tried different extrapolations 
of our results for the 2-body energies and compared to the exact solution. The 



leading corrections from the finite cutoff N scale as 1/ JN + 3/2. A linear 



extrapolation of our results inl/yW + 3/2 works very well. If only the points 
with N > 80 are taken into account, linear and polynomial extrapolations give 
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errors well below 1%. 



2.2 Three and more bodies 



In order to extend this procedure to the three and higher-body sector, two 
additional tools are needed: the Talmi-Moshinsky- Transformation (TMT) [33] 
and Wigner 6j symbols [34]. In the 3-body sector the finite model space is the 
linear hull of the tensor products ni / imi (si) ®4>n 2 i2m 2 (^2) with the restriction 

2 (m + 712) + h + h — N. There are three 2-body contact interactions Vi 

f.i\ 

[i = 1,2,3) and a 3-body contact interaction PQj t . . The contributions of these 



cont • 

interactions are most conveniently calculated in Jacobi coordinates and 
(i = 1, 2, 3). There are three different sets which are labelled by the coordinate 
i not appearing in the definition of . The set s^ 3 , is given by 



^3) 



6 2 



-j= (f 1 - X 2 ) , 
1 _ 



and the other two combinations can be obtained by cyclicly permuting the 
indices. In order to transform the matrix elements between different sets of 
coordinates the TMT is used. The Talmi- Transformation transforms a set of 
coordinates (p, A) to a set (//, A') via 




where d > is a real number. In our case, we need the transformation from 
one set of Jacobi coordinates to another set of Jacobi coordinates. The Talmi 
transformation with d = 3 transforms the set (s^i \ s^) to the set (— s^p, — s^) 
and the set (si 3 \ — §2 ) to the set (— §2 ), respectively. The corresponding 
transformation of the coupled oscillator functions from the first set to the other 
sets is finite, since the total oscillator energy, the total angular momentum, 
and the parity are conserved: 



M L 



E (n' pl l' pl ,n' y l'y]L\n p l p ,n x l x 



n',l',n'l' 

p' p' A' A' 



4>n'^ {?) ® <t>n'lL, (A' 



M L 



The expansion coefficients (n'^l'p,, n x ,l' x ,\ L\n p l p , n\l\^ are the Brody-Moshinsky- 
Brackets. 
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The coupling constants of the 2-body interactions are determined by Eq. (j3J) 
for a given 2-body energy. The coupling constant v ^ of v}^t is determined 
such that a given 3-body energy is reproduced. Using the separability of 
vjolt one finds for the running coupling constant v^ 3 \N): 



ftjUJ7T 3 



I £ ^0,^0,0)/^/^ 

k 



D (k) - 



(5) 



with the eigenvector 

\k)= E z(k;n', 1 l' gl ,n' g / g2 ,L 



N 



nl IL ,n'- Z'_ ,L 

s l s l s 2 s 2 



corresponding to the eigenvalue problem -D(fc) = ( H ° s>c+V ^ +V2+V ' i \\k). 



(6) 



Analogous to the 3-body sector, the model space for four bodies is the linear 
hull of tensor products of HOF with 2 (ni + n 2 + n 3 ) + li + l 2 + l 3 < N. In 
order to determine all contributions of the interactions in one set of Jacobi- 
coordinates, the TMT is used again. The transformation is defined for coupled 
wave functions 



For the TMT related to the mapping {5*2, S3} — > {s' 2 , s'3}, the wave functions 
must at first be recoupled to 



V 
M' 



</>mh ® [0n 2 / 2 (s 2 ) g> (pn 3 h (s 3 )] i23 
The overlap between these wave functions is given by the 6j symbols [31]. 



1-^12 



^13(3 



1 L 




M 





MM' 



-1 



W2 ® <Pn 3 l 3 \ 



1-^23 



(2Li 2 + 1) (2L 23 + 1) 



h h L12 
I3 L L 23 



This construction of the model space can be generalised to A bodies in a 
straightforward way. 
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Figure 3. Energy spectrum for mixed antisymmetric states of three particles in a 
harmonic trap with L p = 1~ as a function of N and exact analytical results from 
[22] . The solid lines indicate the linear extrapolation to N = oo while the dashed 
line indicates the value of the cutoff N at which N ~ §E/{tvjj). 

3 Results for the 3-body sector 



3.1 Unitary limit 



We are now in the position to apply this effective theory to the 3- and 4-body 
sector. First, we focus on the 3-body sector in the unitary limit in order to 
test our method against the exact analytical solution |22j. The corresponding 
ground state energy for the 2-body system is = 1/2 [21J. Since the three 

(2) 

coupling constants v\ are identical, the Hamiltonian is invariant under per- 
mutation of the bodies. It follows that the Hamilton matrix is block diagonal 
in different symmetries of the wave functions. Typical values of the cutoff N 
reached in our study for the 3-body problem are N ph 40. ..70. 

(3) 

For mixed antisymmetric states K ont does not contribute. As an example, 
Fig. [3] presents the spectrum of states with angular momentum and parity 
L p = 1~ depending on the cutoff iV in the unitary limit. Again, the exact 
results are added at N — > oo. The dashed line indicates the value of the cutoff 
iV at which N ph 5E/(hu). It is clear that for the higher excited states an 
extrapolation in N is required. The spectrum can be extrapolated linearly in 
l/\/N + 3 as indicated by the solid lines. 

For completely symmetric wave functions, however, the 3-body interaction 

(3) 

Vcont does contribute. Such wave functions can occur for systems of three 
identical bosons or three distinguishable particles which both display the Efi- 
mov effect. In Fig. HJ we show the spectrum for positive parity and L = as 
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Figure 4. Energy spectrum for completely symmetric states of three particles in 
a harmonic trap with L p = + as a function of N and exact analytical results 
from [22J. Efimov-like states are indicated by squares, other states are indicated by 
crosses. The dashed line indicates the value of the cutoff N at which N ~ 5E/(huj). 

a function of the cutoff N. The 3-body interaction was adjusted such that the 
3-body ground state has energy = e^hjj = —ftw. As shown by Werner 
and Castin [22], there are two different types of states. On the one hand, there 
are states independent of (crosses). On the other hand, there are states 
which depend on (squares) and are called Efimov-like. These Efimov-like 
states are the analog of Efimov states in free space. The exact results [22] 
are given at N — > oo. The dashed line indicates the value of the cutoff N at 
which N rs 5E/(froj). Again, for the higher excited states an extrapolation is 
required. For the non-Efimov-like states a linear extrapolation is appropriate. 
The Efimov-like states, however, show a curvature. In this quadratic 
term has to be included in the extrapolation. Typical extrapolation errors for 
Efimov-like states are of order 2-3% and less than 1% for the other states. 



3.2 Application to e Li atoms 



As an application, we consider a system consisting of three fermionic 6 Li 
atoms in a HOP with a low external magnetic field. Our calculations apply 
to experiments with a few fermions per site in optical lattices. The scatter- 
ing lengths for the three lowest hyperfine states labelled by their hyperfine 
quantum numbers |/,m/),or by integers: |1) = |l/2, +1/2), |2) = |l/2, — 1/2) 
and |3) = 1 3/2 , — 3/2) have been calculated by Julienne [35]. The Feshbach 
resonances occur for magnetic fields near 834 G, 811 G and 690 G. The corre- 
sponding renormalisation energies are determined by the relation between the 
ground state energy and the scattering length a in the 2 particle sector 
(See FigJT]). The range of interactions for ultracold atoms is the van der Waals 
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Figure 5. The energies of first two + states in the 3-body system |1) ® |1) ® |2) 
indicated by circles and pluses (left scale) as function of the magnetic field. The 
solid and dashed lines give the scattering length a\2 from [35J and the van der 
Waals length / V( jw (right scale). 

length i v dW) which is approximately 62.5 for 6 Li, where a# = 0.529 A is 
the Bohr radius. Thus the effective theory is only reliable for regions, where 
a Z 2/ v dw- A discussion of Efimov physics in 6 Li atoms with the three spin 
states |1), |2), and |3) without a trap was given in Ref. |36j. As an example, we 
calculate the 3-body spectrum for the system |1) <g> |1) ® |2) which contains two 
like atoms in the trap. Because of the Pauli principle, there are no Efimov-like 
states and V^t does not contribute. Figure |5] shows the extrapolated energies 
of the first two states with L p = + as a function of the magnetic field (circles 
and pluses). The oscillator length is (3 ~ 65 ag- The solid and dashed lines 
give the scattering length a i2 from [35] and the van der Waals length Z vdW 
(right scale). The hatched areas indicate where the scattering length can not 
be considered large and our theory is not applicable. Outside of the hatched 
areas, our calculation should be accurate up to corrections of order aw/l^w- 
The magnetic field dependence of both states is very weak despite the strong 
variation of the scattering length au- The energies of the ground and first 
excited state are of order 3.2 froj and 5.2 hcu, respectively. 



4 Results for the 4-body sector 

We now turn to the 4-body system. We consider two types of 4-body systems 
both of which can display Efimov physics and focus on the systematics of the 
bound state spectrum. First, we consider two identical fermions with two other 
distinguishable particles. Because of the two identical fermions involved, this 
system is calculationally simpler than the 4-boson system which we calculate 
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Figure 6. Energy spectrum for two identical fermions and two distinguishable parti- 
cles with L p = + in a harmonic trap as a function of N. The system is in the unitary 
limit and = —0.5. Efimov-like states are indicated by circles while other states 
are given by crosses. The solid lines indicate the linear extrapolation to N = oo 
while the dashed line marks the value of the cutoff N at which N ~ 5E/(hw). 

subsequently. Typical values of the cutoff iV reached in our study for the 
4-body problem are N w 20... 30. 

4-1 Two identical fermions and two distinguishable particles 

We start with the case of two identical fermions with two other distinguish- 
able particles. Thus, the wave function is antisymmetric under permutation of 
the identical fermions. Altogether, there are five 2-particle contact interactions 
and two 3-particle contact interactions. The contributions of some interactions 
are pairwise identical because two of the particles are identical. Analogous to 
the 3-body sector, the four coupling constants of the 2-particle interactions 
are renormalised by Eq. (j4]) and the two coupling constants of the 3-particle 
interactions by Eq. (jSJ) for given ground state energies and e^ 3 \ As in the 
2- and 3-body sector, the states, which are independent of e^ 3 ^ can be extrap- 
olated linearly while the Efimov-like states require a quadratic extrapolation. 
Here we require that no minimum exists for 1/N > in order to get stable 
results. In Fig. [61 we show our results for the spectrum states with L p = + 
in the unitary limit and for e^ 3 -* = —0.5. The dashed line indicates the value 
of the cutoff N at which iV w 5E/(hu). For all but the lowest two states the 
condition iV > 5E/(frw) is not satisfied for the values of N reached in this cal- 
culation. As a consequence, a stable extrapolation for N — > 00 is required to 
obtain reliable results. Since the exact results are not known, we estimate the 
uncertainty from the extrapolation conservatively as being equal to the energy 
shift from the last calculated value to the extrapolated value. This procedure 
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Figure 7. Extrapolated spectra for Efimov-like + states in the 4-body sector of 
two identical fermions with two other distinguishable particles (circles) and in the 
3-body sector (squares) for various e^ 3 ) in the unitary limit. The dotted lines give 
our upper bound on the extrapolation error for the two lowest states. The dashed 
lines are guides to the eye. 



should give an upper bound on the extrapolation uncertainty. 



In order to understand the systematics of the trapped spectra, we investigate 
the 4-body spectrum as a function of the input parameters e*- 2 -* and e^ 3 ** analo- 
gous to the Efimov plot in free space [3] . In Fig. [7J the extrapolated spectra for 
two identical fermions with two other distinguishable particles in the unitary 
limit are shown for various ground state energies by the circles. Addition- 
ally, the exact results for the Efimov-like states in the 3-body sector are shown 
(squares). Both the 3- and the 4-body states depend on approximately lin- 
early in this region but the slope of the 4-body states is typically larger. The 
dotted lines give our upper bound on the extrapolation error for the two low- 
est 4-body states. For the higher states the extrapolation error is similar. The 
higher excited states are connected by dashed lines to guide the eye. Note that 
the 4-body ground state is above the 3-body ground state in this case. Around 
e^ 3 -* ~ 0.25, the trajectories of lowest two and of the 4th 4-body state cross the 
next higher trimer state. In free space, the corresponding 4-body state would 
become unstable to decay into the trimer and another particle, but in the trap 
this behavior can in principle be observed experimentally. The other 4-body 
states are different in nature and move parallel with the nearest 3-body state 
and never cross. 
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Figure 8. Spectrum for symmetric + states of four identical bosons in the unitary 
limit for e^ 3 ) = —1 as a function of N. The solid lines indicate the extrapolation 
to N = oo while the dashed line indicates the value of the cutoff N at which 
N « 5E/(hu). 

4-2 Four identical bosons 

We now turn to the system of four identical bosons with L p = + . This system 
is of high experimental interest and the behavior of the states in free space 
is well known [TfffTT] . Figure [8] shows the calculated spectrum in the unitary 
limit as a function of N for = — 1. The dashed line indicates the value of 
the cutoff N at which N ~ 5E/(hu). For most of the states, the condition iV > 
5E/(hu) is not satisfied and we therefore have to rely on the extrapolation. 
Efimov-like states are extrapolated with a quadratic polynomial. Some of the 
extrapolations are shown by the solid lines in Fig. [HJ Again, we estimate the 
uncertainty from the extrapolation conservatively as being equal to the energy 
shift from the last calculated value to the extrapolated value. 

We are now in the position to study the structure of the 3- and 4-body spectra. 
In the original Efimov plot, the 3-body spectrum is studied for fixed 3-body 
interaction while the 2-body energy is varied [Xf3] . Since there is no 4-body 
interaction at leading order, this plot can be extended to the 4-body system 
and has been studied extensively in free space pT6]fTT| . We will compare our 
spectra with the free space results. In Fig. [91 the extrapolated spectra of the 
symmetric + 4-body states for various are shown by the circles. The 3- 
body interaction is fixed by the requirement, that the 3-body ground state lies 
at = — 1 in the unitary limit. Additionally, the 3-body Efimov-like states 
are shown as squares. As before, the dotted lines give our upper bound on 
the extrapolation error for the two lowest 4-body states. The higher excited 
states are connected by dashed lines to guide the eye. Their extrapolation 
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Figure 9. Extrapolated spectra of the symmetric + 4-body states (circles) and 
3-body Efimov-like states (squares) for various e^ 2 \ The 3-body interaction is fixed 
such that e^ 3 ) = — 1 in the unitary limit. The dotted lines give our upper bound on 
the extrapolation error for the two lowest states. The dashed lines are guides to the 
eye. 

error is similar but not shown explicitly. The harmonic confinement has a 
strong effect on the spectrum. Compared to free space, it is no longer true 
that two 4-body states are attached to each trimer state. Moreover, the levels 
appear to interact strongly. There are various avoided crossings of 4-body 
states, e.g. between the 4th and 5th state around e^ 2 -* ~ and possibly also 
between the second and third state. These avoided crossings could be studied 
experimentally by varying e^ 2 ' using Feshbach resonances. 

The dependence on e*- 2 -* can be translated into a dependence on the scattering 
length a using Fig. [TJ For between minus two and zero, the scattering 
length is essentially zero. When e*- 2 -* is varied from zero to one, however, the 
scattering length grows to become infinite at e*- 2 ) = 1/2, jumps to minus infin- 
ity and approaches a negative value close to zero at = 1. This is the most 
interesting region concerning universality and corresponds to the usual Efimov 
plot in free space. In this region, the scattering length is much larger than all 
other length scales and our effective theory is expected to describe systems of 
real atoms with van der Waals interactions. The discrete scale invariance of 
the 3- and 4-body spectra in free space has disappeared in Fig. [9j It would 
be interesting to approach the free space limit by making the trap wider and 
wider in experiment in order to see how the discrete scaling symmetry is re- 
stored. In the theoretical calculation, taking this limit is computationally very 
expensive since the absolute value of the energy cutoff for fixed N vanishes as 
f3 — > oo. Cold atom experiments could serve as a quantum simulator to study 
this issue. 

In Fig. [TOl we show a different variant of this plot. The extrapolated spectra for 
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Figure 10. Extrapolated spectra of the symmetric + 4-body states (circles) and 

3- body Efimov-like states (squares) for various e^ 2 \ For each e^ 2 ) the 3-body inter- 
action is renormalised with the requirement, that the 3-body ground state lies at 

= — 1. The dotted lines give our upper bound on the extrapolation error for the 
two lowest states. The dashed lines are guides to the eye. 

4- body + states are given for various 2-body energies (circles). The 3-body 
interaction is renormalised with the requirement, that the 3-body ground state 
energy lies at e*- 3 ^ = — 1 for each chosen e^ 2 '. Further, the Efimov-like 3-body 
states are added as squares. The extrapolation error estimates for the lowest 
two states are again given by the dotted lines. There are also avoided crossings 
for the higher states near e^ 2 ' ~ 0. In general, the spectrum is very different 
since the 3-body ground state is forced to be at = —1. Experimentally, 
this situation is currently not accessible. It would require a second "3-body" 
Feshbach resonance that can be used to keep the 3-body ground state energy 
constant. 



5 Summary and Conclusions 

In this work, we have studied the universal properties of few-body systems with 
large scattering length in a harmonic trap. We have used an effective theory 
for short-range interactions and identified the cutoff on the harmonic oscillator 
basis functions with the cutoff of the effective theory [31] . The effective theory 
is then formulated directly in the model space and the low-energy constants are 
fixed by matching to few-body observables in the trap. The few-body states in 
the trap are then obtained by diagonalization of the finite Hamiltonian matrix. 
After renormalization, the spectrum is independent of the cutoff up to small 
corrections. If the cutoff can not be taken sufficiently large, it is possible to 
extrapolate to larger values. 
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Our work has focused on systems where the Efimov effect occurs. In this case, £L 
3-body interaction already enters at leading order and the corresponding low- 
energy constant has to be fixed by a 3-body observable. We have checked our 
numerical method against exact analytical results for the 3-body system in the 
unitary limit and found good agreement. Using the magnetic field dependence 
of the 6 Li scattering lengths calculated by Julienne [35], we have predicted the 
lowest two states in a two-state mixture as an example. 

The main part of the work focused on the 4-body system in a harmonic trap. 
We considered two types of 4-body systems both of which can display Efimov 
physics. First, we considered two identical fermions with two other distinguish- 
able particles. Second, we calculated the spectrum for four identical bosons. 
We found a number of interesting features such as avoided crossings of 4-body 
states and 4-body states crossing 3-body states if the 2- or 3-body energies are 
increased. It would be very interesting to observe these features experimen- 
tally. The trap introduces the oscillator length /3 as a new scale and breaks 
the discrete scaling symmetry from free space. In the calculated spectra no 
remnants of the discrete scaling symmetry were observed. 

The extrapolation error of the calculated energies has been estimated conser- 
vatively from the difference of the last calculated value and the extrapolated 
value of the energy. These errors could be improved by going to larger values 
of JV.0 

However, it would also be desirable to understand the leading correc- 
tions from the finite model space in more detail. Similar studies in free space 
have been carried out [37] . 

In general, the trap has a strong influence on the overall structure of the 
spectrum. In particular, the discrete scaling symmetry and the connection of 
two 4-body states to each trimer state known from free space could not be 
observed. In the future it would be interesting to understand this modification 
in more detail by varying the size of the trap. Moreover, the extension to 
more bodies would allow to study the transition from few- to many-body 
systems. While this is in principle straightforward, the computational effort 
grows substantially. 
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